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Abstract 

The beta-decay of a neutron into a bound (pe~) state and an antineutrino in the presence of 
a strong uniform magnetic field (B > 10 13 G) is considered. The beta-decay process is treated 
within the framework of the standard model of weak interactions. A Bethe-Salpeter formalism 
is employed for description of the bound (pe~) system in a strong magnetic field. For the field 
strengths 10 13 G< B < 10 18 G the estimate for the ratio of the bound-state decay rate wi and the 
usual (continuum-state) decay rate w c is derived. It is found that in such strong magnetic fields 
Wb/ W c ~ O.l-j-0.4. This is in contrast to the field-free case, where Wb/w c ~ 4.2 x 10~ 6 [J. N. Bahcall, 
Phys. Rev. 124, 495 (1961); L. L. Nemenov, Sov. J. Nucl. Phys. 15, 582 (1972); X. Song, J. 
Phys. G: Nucl. Phys. 13, 1023 (1987)]. The dependence of the ratio Wb/w c on the magnetic 
field strength B exhibits a logarithmic-like behavior. The obtained results can be important for 
applications in astrophysics and cosmology. 
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I. INTRODUCTION 



It is well known that besides the main decay mode of a free neutron in vacuum 

n — > p + e~ + v e (1.1) 
there is also the neutron decay into a bound (pe~) state (hydrogen atom) and an antineutrino 

n^(pe~)+u e . (1.2) 

Several theoretical estimates for the rate of the latter process have been performed (see 
Refs. 0, 0, [J and references therein). Regardless of the framework employed for the de- 
scription of the hydrogen atom (for example, such as the Schrodinger equation [lj], Dirac 
theory or relativistic Bethe-Salpeter formalism 3]) all the treatments yield the following 
estimate for the ratio of bound-state and continuum-state decay rates: Wb/w c ~ 4.2 x 1CT 6 . 
Therefore, in general one might expect the effect of the neutron bound-state decay (jl.2|) 
to be subdominant if compared with the continuum-state decay (jl.lj) . However, as shown 
below this conclusion does not hold true in the situation where a neutron decays in the 
presence of a strong magnetic field. 

In this work we estimate the ratio Wb/w c in the presence of a magnetic field with strength 
B > 10 13 G. Our study is motivated by a widely accepted view that in diverse astrophysical 
and cosmological environments the physics of neutrinos in strong magnetic fields plays an 
important role. The existence of very strong magnetic fields in proto-neutron stars and 
pulsars is well established. The surface magnetic fields of many supernovae and neutron 
stars are of the order of B ~ 10 12 -r- 10 14 G 
are perhaps as large as B ~ 10 15 -r- 10 16 G 



, |^ . The surface magnetic fields of magnetars 
Very strong magnetic fields are also supposed 
to have existed in the early Universe (see Ref. Q for a recent review). 

The effect of a constant magnetic field on the process (jl.lj) is well documented. In short, 
the constant magnetic field affects the motion of a charged particle in such a way that 
the energy associated with the transverse motion (with respect to the field direction) is 
quantized into Landau levels while the longitudinal motion remains free. The larger the 
field intensity B, the larger the energy separating the Landau levels. Thus, the decay rate 
w c exhibits the following dependence on the field intensity B. For < B <C B cr , where 
B cr = (A 2 — m 2 )/2e = 1.2 x 10 14 G with A = m n — m p being the nucleon mass defect 
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(hereafter we use the units h = c = 1), the effect of a magnetic field on both an electron 
and a proton is small and w c remains practically insensitive to B. As the field intensity 
approaches the value B cr (note that B cr > B e , where B e = m 2 e /e = 4.414 x 10 13 G is the so- 
called Schwinger field) the effect of Landau quantization on the transverse electron motion 
becomes considerable u c ~ B 8]. If B cr < B B' cr , where B' cr = [(m n — m e ) 2 — mi]/2e = 
1.25 x 10 17 G, the electron (due to the energy conservation law) can occupy only the lowest 
Landau level, and it was shown in e 8] that the decay rate grows as w c ~ B. In the case of a 
superstrong magnetic field (B > B' cr ) the proton can occupy only the lowest Landau level 
and one has again a monotonic dependence w c ~ B j2| . Finally, as the field strength exceeds 
the value B ~ 1.5 x 10 18 G the modification of the strong forces, which bind the quarks in 
nucleons, is such that it can close the mass gap between the neutron and proton making the 
neutron stable Q (see also Q). In what follows we concentrate on field strengths in the 
range B e < B <C B p , where B p = rrv^/e = 1.5 x 10 20 G is the Schwinger field for the proton. 

To our knowledge, no theoretical analysis has been published for the process (jl.2)l in a 
(strong) magnetic field. It should be noted that such study is hampered by the fact that 
in calculating the decay rate Wb one must know the eigenenergies and eigenstates of the 
hydrogen atom in an external magnetic field. In the case of a strong field the latter problem 
has been studied in a number of papers, where numerical and semianalytical results have 
been derived within different approaches based on the Schrodinger 
H Q H Dirac HQQ, and Breit 3| 

equations. A common result of all the 



treatments is that, compared to the field-free case, the ground state of the hydrogen atom is 
(1) very tightly bound and (2) well separated in energy from excited states. The shape of the 
atom is also affected by the field. The transverse atomic size is approximately determined 
by the magnetic length <xb = 1/y/eB <C ao, where ao = (e 2 m e )' 1 is the Bohr radius, and the 
longitudinal atomic size an depends on the binding atomic energy. In particular, in the case 
B > B cr one has for the ground state a cigar-like shape with an <C a , while for the excited 
states a\\ ~ a j^. It should be noted that these remarkable effects of a strong magnetic 
field on the hydrogen atom are obtained assuming the absence of transverse motion of the 
atom as a whole. However, even at a moderate transverse atomic velocity a rather strong 
electric field is induced in the rest frame of the atom pushing the electron and the proton 
apart. As a result, the binding energy decreases as the transverse velocity increases. Thus, 
in contrast to the field-free case, the center-of-mass atomic motion does not decouple from 
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the relative internal motion, and therefore a correct description of the hydrogen atom in a 
strong magnetic field should involve the two-body approach jssl . I^j ] . 

By analogy with the field-free case jsj, the ratio Wb/w c can be roughly estimated using a 
phase-space argument that does not depend on a formal theory of weak interactions. The 
phase-space volume available for the final products in the process (II. lj) can be presented 
by the function f(B, A) that depends on the field strength and nucleon mass defect. For 
the bound-state decay (jl.2|) . the available phase-space volume is given by £^|-?/>(0)| 2 , where 
e v is the neutrino's energy and i/j(0) is the ground-state wave function of a hydrogen atom 
evaluated at zero distance between an electron and a proton. The ratio of the bound-state 
and continuum-state decay rates is approximately equal to the ratio of the corresponding 
phase-space volumes 

For example, using the phase-space argument in the field-free case (B = 0) we get the result 

W b .6 ( A t \ _on„ in-6 



~ vreM 1 » 2.9 x 10 _b , (1.4) 

which is of the same order of magnitude as the accurate calculations of Refs. |l|, |2j, |3j). In 
the case of a strong field (B cr < B <^ B' cr ) we have 

f(B, A) ~ eBml/ir 2 and |^(0)| 2 ~ (vra^l)- 1 = (Tra^/eB)- 1 . 

Using the estimates for an reported in Ref. j^j, we might expect [in accordance with Eq. (|1.3j) ] 



that 



2 /A \ 2 

Wfe 7re ao / A x 



1 >0.1. (1.5) 

w c ay \w^ e / 

This result is almost five orders of magnitude larger than in the field-free case (jl.4j) and thus 
it calls for a more detailed and rigorous theoretical analysis, such as carried out below. 

The paper is organized as follows. In Sec. |H]a general formulation of the problem is given 
in the context of the standard quantum field theory employing a Bethe-Salpeter formalism. 
Then, in Sec. IIII1 we discuss in detail a structure of the wave function of a bound (pe~) state 
in a strong magnetic field. Further, specific approximations to the Bethe-Salpeter equation 
are developed in the cases B cr < B -C B' cr and B' cr < B <^ B p . Sec. HXlis devoted to the 
derivation of an estimate for the bound-state decay rate. In addition, the asymptotic formula 
for the ratio of bound-state and continuum-state decay rates is obtained. The conclusions 
are drawn in Sec. M 



II. GENERAL THEORY 



In the case of field strengths relevant to this work {B e < B <C B p ) we can neglect the 
effect of a magnetic field on both the propagator of W boson and the weak form factors. 
Therefore, the transition matrix element for the decay process (jl.2|) can be written as follows 



T fi = — I x E (x, x) 7m (1 + a l5 )ij n ^(l + 75 )<M 4 :r, (2.1) 



G 

where ip n and ip u are the Dirac wave functions of the neutron and neutrino, respectively, 
G = Gfcos6 c , 6 c is the Cabibbo angle, and a = 1.26 is the ratio of the axial and vector 
constants. The conjugate \ E °f the wave function \ E describing the bound {per) system 
with energy E is determined as \ E = 7°X 7°- m the ladder approximation the wave 
function \ E satisfies the following Bethe-Salpeter equation 

X E {x e ,x p ) =ie 2 J G e {x e ,x' e )^G p {x p , x' p )Y D, u {x' e , x' p ) X E {x' e , x' p )d%d% (2.2) 

Here D,„,( the photon propagator. The electron and proton propagators satisfy the 

Dirac equation in an external magnetic field 

[YHd^T eA,) + m e/p }G e / p {x,x') = 5^{x - x'). (2.3) 

In what follows we use the symmetric axial gauge of the vector potential 

A(r) = ~[Bxr]. (2.4) 

The photon propagator is supposed to be invariant under the time and space translations 

D^{x,x') = D pu {v-v',t-t') } (2.5) 

whereas the electron and proton propagators are invariant under the time translation 

G e/p {x, x') = G e/p {r, r', t-t'). (2.6) 

Their properties with respect to the space translation are determined according to Eqs. ()2.3|) 
and ()2.4|) as follows 

G e/ P{Y _ rC) Y > - rcj t- t') = exp {±^[B x r c ](r - r')} G e / p (r, r', t-t'). (2.7) 
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From (J23J), (|2~7I) . and (l2~2l we receive 

X E (x e , x p ) = exp(-iET) exp jz (p + |[B x r]) r} r^ p (r, t), (2.8) 

where P is the total pseudomomentum which corresponds to the total momentum in the 
field- free case , t = t e — t p and r = r e — r p are respectively the relative internal time and 
space position, while T and R refer to the center-of-mass time and space coordinates. 
Thus, passing in (|2.2jl to Fourier-transform according to the rule 

G^(r,r',t-t') = J G^(r,v',e)exp[-ie(t-t')}^, (2.9) 

we obtain the equation 

v E > F (r,t) = ^ / exp[-ie(t - t')}G e (r, R + r', e)^G p (0, R, E - e)^ v D ^{r' , t') 
lit J ~ 

x exp {i (P + |[B x r']) R'} r] E ' p {r' ,t')dedr'dK'dt' . 

(2.10) 

We will assume that the vector potential (|2.4|) is specified in the rest frame of the neutron, 
which is supposed to be not affected by a magnetic field. Thus, the neutron wave function 
is given by 

^ n = u n e- m -\ u n =^j (u;>„ = l), (2.11) 

where p n = (m n ,0), w n is the two-component neutron spinor. The neutrino is supposed to 
be massless and not interacting with a magnetic field. Therefore its wave function can be 
presented as follows 



^ = ^e-^- p " r \ u„ = 4?l ! <'V-'V : 1 ) ■ (2.12) 



V2 



where p u = (e u ,p v ), e u = \p u \, and v u is the two-component eigenspinor of the spin- 
projection operator, (o-p u )v v = — e v v v . 

Accounting for Eqs. ()2.8|) . (j2.11|) . and (|2.12j) . the matrix element (|2.1|) takes the form 

T fi = -^(2vr) 4 5(m n - e„ - E)8®{p v - P)rf ' p (0, 0)^(1 + a 75 K 7 "(l + 75 )u„. (2.13) 
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Using the rules 1 \2n5(E)\ 2 = 2n5(E) and |(2tt) 3 5( 3 )(p)| 2 = (2n) 3 5^(p), we get for the rate 
Wb of the decay process (jl.2j) 

r2 



W b 

where 



M = rf- F (0, 0)7„(1 + a7s)ti„7"(l + 7sK (2.16) 



is a reduced matrix element, the sum s ^} pe ) runs over all bound states of the (pe ) system 
and the average over neutron spin states is assumed. 

III. THE BETHE-SALPETER WAVE FUNCTION 

The key element that determines the decay rate (j2.14|) is the wave function of the bound 
(pe~) system. We solve the Bethe-Salpeter equation (j2.1U|) neglecting the retardation ef- 
fects |22| (this amounts to the non- relativity of internal motion), so that the photon propa- 
gator in the Coulomb gauge assumes the form 

D^{T,t) = (3.1) 

Inserting (|3.1j) in the Bethe-Salpeter equation ()2.1()|) . we obtain 

^' P (r, 0) = [ G e (r, R' + r', eh°G*>(0, R', E - e)^-, 

In J r 

x exp ji (P + |[B x r']) R'} r/ £ ' p (r', 0)dedr'dR'. (3.2) 

To carry out an integration over e in the kernel of this equation, we expand the electron 
and proton propagators into series over Dirac eigenstates in a magnetic field keeping only 
the positive energy pole contributions [28]: 

G « /P(rir , eHE «W^|M (3 . 3) 



If z axis is directed along the magnetic field vector B, then e K = y m e/p + V 2 Z + 2neB is the 
energy of the Dirac eigenstate specified by a set of the quantum numbers k = {n,p z ,j z , s}, 
where the discrete numbers n = 0, 1, 2, . . . denote the Landau levels, p z is the longitudinal 



1 The normalization length and time are set to unity. 
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momentum, j z is the z component of the total angular momentum, and s = ±1 is equivalent 
to the spin quantum number in the field- free case (see, for instance, 2^|). In cylindrical 
space coordinates, the Dirac wave functions in Eq. (|3.3j) can be presented as the products 
of longitudinal and transverse parts 



(*), 



e/p\ 



(3.4) 



where X*)(p z ,p) is the transverse spinor labeled by a set of the quantum numbers k 



{n,j z , s}. The explicit forms of the transverse electron and proton spinors are 

/ y/m e + £« e 0„_ij z -i/2(p) \ 



Xf^ +1 \p z ,p) 



Pz 



>n-l,j z -l/2 



2neB 



/ o \ 

Vm e + £ Ke ^ + i /2 (p) 



2neB 
m e + e K( 



0n Jz +l/ 2 (p) 



X^> +1 \ Pz ,p) 



( \/ m P + £ K p <f>n,j z -l/2(p) \ 



Pz 



y/ m p + e Kp 
2neB 



>n,jz-l/2 



(p) 



m p + e Kp 



X^'^ip^p) 



n-l,j z +l/2(p) 

I o N 

^/m p + e Kp (j) n - 1Jz+1 / 2 (p) 



2neB 



m p + e Kp 

-Pz 



where the transverse Landau orbitals are given by |2l|, |2 

'n p ,m\P 



?n-l Jz +l/2(P) 



/ 



2ira 2 B (n p — m)\ 



exp(im(p) exp 



P 



4a 2 B J \2a B 



n 2 \~ m/2 I rr 

P \ r(-m) / _P_ 



(3.5) 



2a% 



(n p = 0,l,2,..., 



n p - m = 0,1,2,...) (3.6) 



with Ln p m ^ being a generalized Lagguerre polynomial. 
Using Eqs. ()3.3|) and (|3.4|) . we have 



^ , J e- y/e% +p 2 z + i0 



where = ^Jm 2 e j p + 2neB is the energy of the transverse motion. Inserting ()3.7|) in ()3.2j) 
and integrating over e, we receive 

1 E-Jel + q l-Jel+{P z -q : V 



x 



\f*>{P z - g „R' ±)7 =L=exp {i ( P± + |[B x p']) Rl} 



xr/ B ' p (r', 0)dg z dp'd/dRi,(3.8) 

where the sum over k p involves only the transverse proton spinors with j z = ±1/2 [according 
to Eqs. (|3.5J) and ()3.6|) the transverse proton spinors with j z ^ ±1/2 vanish at p = 0]. 

We consider below the following two cases of field strengths: (1) B cr < B <^ B' cr and 
(2) B' cr < B <^ B p . For each of these cases we develop the corresponding approximation to 
the Bethe-Salpeter equation (J3.8j) taking into account the kinematical regime of the decay 
process f)1.2jl . 



A. The case B cr < B <C B' cr 

Recall that in the case B > B cr the energy conservation law dictates that the electron 
in the usual decay process (|1.1|) can occupy only the lowest Landau level. Considering 
the bound-state decay process (jl.2j) in the case B > B cr , we note that [as follows from 
Eq. ([3.8p ] the electron can (virtually) occupy not only the lowest Landau level. However, 
the probability for the electron to occupy the excited Landau levels is vanishingly small due 
to the following fact. Virtual electron transitions between the lowest and excited Landau 
levels are induced by the Coulomb electron-proton interaction, which in the present case of 
field strengths is much weaker than the interaction of the electron with a magnetic field. 
Thus, the electron contribution from excited Landau levels to the wave function of the bound 
(pe~) system is negligible. Therefore, one can leave in Eq. ()3.8|) only those electron terms 
that correspond to the lowest Landau level (this amounts to the adiabatic approximation). 
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In accordance with Eq. (|3.5|) . the transverse electron spinors for the lowest Landau level 
(n — 0) are given by 

A? ,i "~ 1) (Pz, p) = u(p z )(t> Qim (p), m = j z + ~ = 0, -1, -2, . . . , 



/ 







2£ Ke 


~Pz 



(3.9) 



where e Ke = y / mf"+~pf. In addition, the non-relativistic approximation applies to the proton 
transverse spinors and energies: 

A?' lA+1) fe, P) = 4 + Vn, (p), ^'~ 1/2 '~ 1} (Pz, P) = 4"Vn-l, (p) [<t>-lM = 0], 



£k p = m p + nu p + u p +) 






W 



(-) 



1 



w 



(3.10) 



where Wp = eB/m p is the Larmor frequency for the proton. 

Taking into account Eqs. ()3.9|) and ()3.10|) and using nj o(O) = y^eB/2n, we arrive at 



T] E ±> p (r,0) = f(r)u^ 



v ' 



(3.11) 



where .EL = E + + uo v and 



fk(r) = ~ 




2tt V 2tt 



^(j>0,m(p 



exp[iq z (z - z')]u k (q z )ui(q z 



%,m(.P' + R±)0n,o(-R± 



v 7 ^ 



L= exp (p ± + ~[B x p']) R' ± } f^dq.dp'dz'dR^. 

(3.12) 



Here 



^e(g«) = \/ml + q 2 z and E^n, - = m p + ntUp + — 



k,l = 1,2,3,4 stand for the spinor indices and a sum over the repeating spinor index I 
is assumed. To carry out an integration over R' ± in the kernel of Eq. (|3.12jh we use an 
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expansion of the undisplaced Landau state over an infinite set of the displaced Landau 
orbitals Q 

K,m{p) = exp («|[B X p]p^j ^ (j) n -m',m-m' (pa)<f> n ,m' (P ~ Po)- (3.13) 

m' 

Involving the following integral 

f 2ir 

/ eX P(* P ±p)0^- m ',m-m'(p)0n-^,o(p)^P= JgK-m' ,N-m>(Pc)<Pn~ m ,N~m{pc) , (3.14) 

where pc* = e[P^ x B]/(e.B) 2 , and introducing 

</(r) = exp /(r - p c ), (3.15) 

we get 

/ \ e2 ST a f„\ f e M i( lz(z- z')K(g^Kfe) 

gi(v')dq z dp'dz' . (3.16) 



m 



y/tf-pcY + z*' 

Let us seek a solution to this equation in a form of the following expansion 

#( r ) = h ( m] z )^o,m(p)- (3-17) 

rn 

Inserting it in Eq. (|3.1fij) we obtain an infinite set of the coupled integral equations 

h ( \ 1 f ex Pfe(^ - z')]u k (q z )ui(q z ) , , 

h h (m;z) = — > / — rV mrn >(pc,z)hi{m;z)dq z dz . (3.18) 

Zir^J E + - E e (q z ) - E p (\m\, P z - q z ) 

Here the one-dimensional potential functions V mm r(pa, z) are given by 

t r / \ 2 /" 0O,m(p)0O,m'(p) 

^'^,^ = -6 / 7= = 9 ^P- (3.19) 

■/ V(p-Pc) +^ 

Examining their properties, we deduce that in the case Pj_ = the infinite set of integral 
equations ()3.18j) becomes decoupled 

h(m; ,) = J- / ^f-^MM vM)hl(m , J )dqM , (3.20, 
2tt J E + - E e (q z ) - E p (\m\, P z - q z ) 

where 



|m| 1 

2 / /v2 \ / ^2 \ 2 4 / —2 



Mz) = -e / rfp = -e exp — -=- — =- W ' 



(3.21) 
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with Wfj,, v being a Whittaker function 

In the non-relativistic limit, the equation (|3.20|) is equivalent to the one-dimensional 
Schrodinger equation for a hydrogen atom in the adiabatic approximation 

1 d 2 



V m (z) + El 



ipm(z) = 0, (3.22) 



2fi dz 2 

where /i = (m e + m p )/2m e m p is the reduced mass and 

Ef n = E + -m p -m e - — — — r - \m\u p , ip m (z) = exp ( lPz7neZ J h 2 (m; z). (3.23) 



2(m e + m p ) 

Here it should be pointed out that the use of the non-relativistic approximation (|3.22)l for 
the treatment of the bound states in the case of a magnetic field with strength B > B cr 
can be justified for the following reasons j^jj: (1) the shape of the transverse electron wave 
function determined by the Landau orbital (|3.6|) in the relativistic theory is the same as 
in the non-relativistic theory; (2) the electron remains non-relativistic in the z direction as 
long as the binding energy \E^ \ <C m e . 



12. 



The one-dimensional Schrodinger equation ()3.22|) has been studied in 
Here we outline its basic properties that are important for our purposes. The eigensolutions 
of Eq. ()3.22|) can be categorized in two distinct classes: (1) the states having no node and (2) 
the states having node(s) in their wave functions. The states having no node in their wave 
functions are tightly bound. For each m there is one such state, which is the most tightly 
bound (the ground state) in the case m = 0. The ground-state energy E" <C —13.6 eV, 
i.e. it is much lower than that of a hydrogen atom in the field-free case. On the contrary, 
the states having nodes in their wave functions are weakly bound. For example, in the case 
m = the state with one node in its wave function has about the same energy as the ground 
state of a hydrogen atom in the field-free case, E" ~ —13.6 eV, whereas the states with 
more than one node in their wave functions have the higher energies (E« > —13.6 eV). 



B. The case B' cr < B <C B p 

When considering the influence of very strong magnetic fields on the neutron beta-decay, 
one should be careful about the effect of a magnetic field on the anomalous magnetic mo- 
ments of the neutron and proton. In particular, it is known (see, for instance, Ref. Q) that 
the interplay between the anomalous magnetic moments of the neutron and proton shifts 
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the energies of these particles making the neutron stable in the case B ~ 1.5 X 10 18 G. How- 
ever, if the field strength is B < 1.5 x 10 18 G the effect of the nucleon anomalous magnetic 
moments is still subdominant. Note that the corresponding shift of the electron energy due 



32, 



(see also 



to the electron anomalous magnetic moment vanishes in such a strong field 
Ref. Q). In what follows, we neglect the effects connected with the anomalous magnetic 
moments, assuming that the field strength does not exceed the value B ~ 1.5 x 10 18 G. 

In the case of a superstrong magnetic field (B' cr < B <C B p ) both the electron and the 
proton in the usual decay process (|1.1|) can occupy only the lowest Landau level. Using the 
same arguments as in the preceding subsection with regard to the electron's dynamics in the 
bound-state decay process (jl.2j) . we again arrive at Eq. (J3.18)) . However, noticing that the 
(virtual) proton transitions between Landau levels induced by the Coulomb electron-proton 
interaction are almost completely suppressed, we can neglect the couplings between different 
values of m in Eq. f)3.18j) and consider only the case m = 0. Thereby, instead of the infinite 
set of coupled equations ()3.18|) . we obtain the single integral equation 



hk{0]z) = — 
Ztt 



exp[iq z (z - z')}u k {q z )ui{q z ) 



E + -E e (q z )-E p (0,P z -q : 
The corresponding Schrodinger equation is 

1 d 2 



■Vo (pc,z')hi(0;z')dq z dz'. (3.24) 



2/i dz 2 

where pc plays a role of parameter and 

p2 



V 00 (pc,z) + El(p c ) 



ipo(pc,z) = 0, 



K{p c )=E, 



m p — m e 



2(m e + m E 



exp 



iP z m e z 



(3.25) 



h 2 {0;z). (3.26) 



Note that if pc -C a#, we have 
V o(pc,z 



VnZ) 



— \ — erfc 
a B y 2 



exp 



2a% 



7r e 
2a B ' 



e 
\z\ 



where erfc(?/) is the complementary error function 



0. 

oo, 
(3.27) 



IV. THE DECAY RATE 



In this section, using the results of Sec. IIHI for the wave function of the bound (pe ) 
system, we derive the bound-state decay rates in the particular cases, namely, B cr < B <^ 
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B' cr and B' cr < B <C B p . Then we utilize the obtained bound-state decay rates in the 
extrapolation procedure, when treating the general case B e < B <C B p . 



A. The case B cr <B«B[ r 

Using Eqs. (|3.11|) . (|3.15|) . and f)3 . 1 7j) we perform an integration over P in Eq. ()2.14|) and 
obtain for the bound-state decay rate 

i2 



w b 



G f 



+6[m n -e v - E_Ap»)]\M-APv)\ 2 } T^p- (4-1) 

Here the indices r label the solutions of Eq. (J3.18|) with corresponding eigenenergies E± tT (p u ) 
[E- tT (p u ) = E +tT (p u ) +Up]. In accordance with Eq. ()2.15j) . the reduced matrix elements are 
given by 

M±, T (p v ) =X)^,mW[«i ±) 7M( 1 + aj 5 )u n ][h T {m; 0) 7 "(1 + 75K], (4.2) 

m 

where p u = e[p u x B}/(eB) 2 . 

To evaluate the reduced matrix element (|4.2|) . we note that for the field strengths under 
consideration one has p™ ax /as = p™ ax aB < 1, where p™ ax is the maximal possible value for 
the neutrino's momentum in the process (jl.2j) . This allows us to put p u = in Eq. (|4.2j) . so 
that we get 

eBG 2 



w b 



E / {*[A-e v -m e - J E|(R,,,)]|>f + ,r(p^)| : 



+5[A - e „ - m e - Wp - 4(p^)]|^_, r ( Pl/ )| 2 } (4.3) 



Here 



M +jr ( Pi/ ) = v^[^: 2 z *(o)-^r*(o)][(i+«)^ n , 1 ^ 2 -2«^ 2 ^ il ], 

M_, T (p,) = >/2[/i%"*(0) - ^;X*(0)](1 - a)u n>2 v u>2 , (4.4) 
where in accordance with Eq. ()3.18|) the wave functions h^' z are the solutions of the equation 

CW4/ expg^-^f^ Vi(,0 ft r-(,')^^.(4.5) 

2tt 7 E"(p u<z ) +m e + m p - E e (q z ) - E p (0,p w , z - q z ) 
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In the non-relativistic approximation 

2 

E\}(p VtZ ) = El + Pv * V ^"(0)=^r(0), ^r(0) = 0, (4.6) 

where ip T are the solutions of the Schrodinger equation (|3.22|) for m = with energies 
Eq t = Et. The index r = 0, 1, ... corresponds to the number of nodes in the wave function 

</v0). 

Using Eq. (|4.6|) and neglecting tu p and p 2 /2(m e + m p ) with respect to A — m e — E", we 
obtain (after averaging over neutron spin states) the following result for the bound-state 
decay rate: 

w b = ^(1 + 3a 2 ) l^(0)| 2 (A - - ^) 2 , (4.7) 

T 

where the sum runs over even states (r = 0, 2, . . .). The odd states (r = 1,3,.. .) do not 
contribute to the bound-state decay rate because their wave functions vanish at the origin 

[Vr(0)=0]. 

Let us derive the ratio of the bound-state decay rate Wb [see Eq. (|4.7j) ] and the continuum- 
state decay rate w c , i.e. the decay rate of the usual process (jl.lj) . To facilitate an accurate 
estimation of the ratio Wf,/w c , we calculate w c under the same assumptions that have been 
made to obtain Eq. (|4.7|) . i.e. we assume p v = and uj p , p 2 /2(m e + m p ) < A - m e . Thus, 

we get 

12 / A _ lA X ' 



where 



(4.8) 



(7 = \ (— 2 + 2) - l _ A arccosh (A\ . (4.9) 

It is useful to compare the result (|4.8jl with that in the field-free case (see also Refs. 

HQ): 

™6 2vr 2 E„l^n(0)| 2 fA-e n 



w c C ml \ m, 



1 , (4.10) 



where 



C 1 f 2A4 8^ 1+ A h( A ) (111) 

60 \ m* m 2 , / y m 2 4m e \m ( 



n = 1,2,... is the principal quantum number, ip n [|?/>n(0)| 2 = {<&m z afy~ x ] and e n = 
— (2aon 2 ) -1 are respectively the wave function and the binding energy of the ns state of 
a hydrogen atom. It is seen that both expressions (|4.8|) and (|4.10j) have a similar structure. 
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However, the important difference consists in that ip T are the one- dimensional wave functions 
while ip n are the three-dimensional ones. This difference is also reflected in the appearance 
of the factor m~ l in the right-hand side of Eq. (J4.8)) instead of the factor m~ 3 occurring in 
the field- free case (J4.1(Jj) . 



B. The case B' cr < B « B p 

By analogy with the previous subsection, in the case of a magnetic field with strength in 
the range B' p < B <C B p we have 

% = y^ J 2n5[m n - e v - E + , T (p u )]\M + Ap,)\ 2 ^ (4.12) 
where the reduced matrix elements are given by 



M+,rM = exp [u p +) l,(l + oc lh )u n \ [h T (0; 0) 7 "(1 + 7sK] (4.13) 

with h T being the solutions of Eq. (|3.24j) . Following the approximate procedure developed 
in the preceding subsection, we deduce that 

eBG 2 



w b 



8tt 2 



1 + 2a + 5a 2 ) £ |^ T (0)| 2 (A - m e - E^) 2 , (4.14) 



where ip T and E} are specified in ()4.6jl . And for the ratio of bound-state and continuum-state 
decay rates we obtain 

2 

(4.15) 



w r Cm P \ m P 



It is remarkable that this expression is identical with the one obtained in the case B cr < 
B <^ B' cr [see Eq. ([4.8)1 ] . We utilize this fact below, when estimating the ratio Wb/w c in the 
general case B e < B B p . 



C. The case B e < B <C B p 

A smooth dependence of w c on the field strength B ^,0] combined with the identity of the 
expressions ()4.8|) and ()4.15|) allows us to extrapolate the obtained results for the ratio Wb/w c 
in the whole range of considered field strengths. Specifically, we assume the expression ()4.8|) 
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1)4.15)1 ] obtained in the case B cr < B <^ B' cr [B' cr < B <^ B p ) to be valid for estimating 
the ratio Wb/w c in the general case B e < B <C B p , which incorporates the above particular 
ranges of field strengths. 

Let us note that in the field-free case (B = 0)_t 
to the ratio Wb/w c [see Eq. ()4.1(J)) ] is about 20% 
to the field-free situation, in the presence of a strong magnetic field the excited states are 
very weakly bound in comparison with the ground state, we neglect their contribution to the 



le relative contribution of excited states 
3|. Taking into account that, in contrast 



ratio Wb/w c . The ground-state wave function ip(z) can be approximated as follows [12L 134] 

4>{z) = (2 /U |E!l|) 1/4 exp f-y/^mWz^ , (4.16) 
where E" is the ground-state energy. Thus we get {ji ~ m e ) 



Wb 

w c 




2 

(4.17) 



In the field-free case the corresponding result is given by (accounting only for the ground- 
state contribution) 

* jMKflzl-!)', (4.18) 

where e = — e 4 m e /2 is the ground-state energy of the hydrogen atom. 

Fig. [U displays the numerical results for the ratio Wb/w c calculated in accordance with 
Eq. ()4.17jl usin g tw o asymptotic formulas for E^: (1) the well-known formula from Ref. 
(see also Refs. Q, Q, Q) 

E W=e\n 2 (-^-\ (4.19) 



(4.20) 



y 2m e \e\ / 

and (2) the recent formula from Ref. |20| 

2 f c 00 e J B/2m e |e| 



E" = eln z 



ln 2 {eB/2m e \e\)_ 

where Cqq = 0.2809. It is seen that both asymptotic formulas give the same order of 
magnitude for the ratio Wb/w c . In both cases the results exhibit a linear dependence on 
\og lQ (B/B e ). However, the results obtained on a basis of Eq. 1)4.19)1 are approximately two 
times larger in magnitude than those obtained on a basis of Eq. ()4.20j) . Due to the fact that 
the asymptotic formula ()4.20j) yields more accurate values for the ground-state energy 
than the formula 1)4.19)1 does (see Ref. j^] for details), we can conclude that the estimate 
for the ratio Wb/w c using Eq. ()4.20)) is more realistic. 
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V. SUMMARY AND CONCLUSIONS 



In summary, we have considered and analyzed theoretically the neutron decay into a 
bound (pe~) state and an antineutrino in the presence of a magnetic field with strength 
B > 10 13 G. The amplitude of the bound-state decay process has been formulated within 
the framework of the standard model of weak interactions. For the description of the bound 
(pe~) system we have employed the Bethe-Salpeter equation. The approximations to the 
Bethe-Salpeter wave function of the bound (pe~) state have been specified depending on the 
field strength B. We have derived the bound-state decay rate Wb in two particular cases, 
namely, B cr < B <^ B' cr and B' cr < B <^ B p . In both cases, the identical expressions for 
the ratio of bound-state and continuum-state decay rates Wb/w c have been obtained. We 
have estimated the ratio Wb/w c in the general case B e < B <C B p using two asymptotic 
formulas j3, for the ground-state energy of the hydrogen atom in a strong magnetic 
field. For both asymptotic formulas a logarithmic-like behaviour Wb/w c ~ a log 10 (B/B e ) + b, 
where a and b are positive constants, has been determined. 

The numerical estimate for the ratio of bound-state and continuum-state decay rates has 
been performed. It has been found that in contrast to the field-free case, where the bound- 
state decay m ode is suppressed by a factor of about 4 x ^ flflQ as cotupared with the 
usual (continuum-state) decay mode, in the presence of a strong magnetic field B > B e the 
ratio Wb/w c is of the order O.l-i-0.4. This remarkable finding can be important for the physics 
of supernovae and neutron stars, where magnetic fields with strength B > 10 13 G may 
exist. In particular, the high value of the neutron bound-state decay rate Wb can influence 
the nucleon balance in protoneutron stars and the hydrogen fraction in the atmosphere of 
neutron stars and magnetars. The high value of Wb can affect the neutrino spectra from 
astrophysical objects with strong magnetic fields, because in the neutron bound-state beta- 
decay process the energy distribution of antineutrino is peaked about A — m e . However, for 
estimating these possible astrophysical effects one should consider a more involved problem, 
namely, the bound-state beta-decay of a neutron moving in the presence of a strong magnetic 
field and dense matter (for instance, such as a neutron star). This implies modifications of 
the photon, electron, and proton propagators in the Bethe-Salpeter equation ()2.2|) due to 
many-particle effects 3|. In addition, an electric field induced in the rest frame of a neutron 
(due to its transverse motion in the matter rest frame) should be accounted for, because 
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for the transverse neutron velocities v± ^> obI^'I the induced electric fields are as strong 
as to pull the hydrogen atom apart. Note, however, that in dense matter an electric field 
can be strongly screened by the surrounding medium. The screening also modifies the 
electron-proton interaction. In the context of the present analysis the latter factor becomes 

— 1/3 

appreciable if ag > l s , where l s is a screening length for the medium. Since l s ~ n c , 
where n c is a density of charged particles (electrons or protons) in matter, we obtain the 
following criterion a# > ric 1 ^ (or B < n'J 3 /e) which indicates the situation where the role 
of screening can not be neglected (such situation can be encountered, for example, in the 
interiors of a protoneutron star). 

Let us remark that the high value of the bound-state decay rate Wf, can be also important 
for cosmological applications. If the controversial hypothesis of strong magnetic fields influ- 
encing the usual beta-decay process in the early Universe |2( is realistic, then one has a right 
to expect that the neutron beta-decay into a bound state of the (pe~) system might have 
substantial consequences for big-bang nucleosynthesis and the production of light elements 
in the early Universe. 

Note that in the terrestrial laboratory environment the strongest magnetic field that can 
be produced is of the order of 10 7 G (see, for example, j^) which is much lower than B e . 
However, the results of our present analysis can be used even in the case of such fields if the 
neutron beta-decay takes place in semiconducting or dielectric media, where a small effective 
mass for the electron and a large dielectric constant reduce the Coulomb force relative to 
the magnetic force (31I ] . 

Finally, it is straightforward to generalize the results obtained in this work to the bound- 
state beta-decay of a nucleus with charge Z — 1 in a strong magnetic field. This is realized 
by replacing the nucleon mass defect A with the corresponding value for the nucleus under 
consideration and taking into account that the atomic energy and radius are given by e = 
—Z 2 e 4 m e and a = (Ze 2 m e )~ l , respectively. After performing such a procedure, it can be 
deduced from Eqs. (pLTTj) . (ETTHj) and/or (jOUJ) that w b /w c ~ Z[alog 10 (Z- 2 B/B e ) + b], where 
the positive constants a and b do not depend on Z. Thus, the dependence of the ratio Wb/w c 
on the nuclear charge is different from that in the absence of a magnetic field, in which case 
Wb/w c oc Z 3 . 
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FIG. 1: The estimate for the ratio of bound-state and continuum-state decay rates in the case of 
a strong magnetic field. The cases, where the asymptotic values for E" are derived according to 
Ref. 35] [see Eq. gUJ)] and Ref. Q [see Eq. ()4.2(Jj) ]. are represented by the solid and dashed 
lines, respectively. 
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